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ABSTRACT 

The Be/X-ray binary A 0535+26 showed a normal (type I) outburst in August 2009. It is the fourth in a series 
of normal outbursts associated with the periastron, but is unusual by presenting a double-peaked light curve. 
The two peaks reached a flux of ~ 450 mCrab in the 15-50keV range. We present results of the timing and 
spectral analysis of INTEGRAL, RXTE, and Suzaku observations of the outburst. The energy dependent pulse 
profiles and their evolution during the outburst are studied. No significant differences with respect to other 
normal outbursts are observed. The centroid energy of the fundamental cyclotron line shows no significant 
variation during the outburst. A spectral hardening with increasing luminosity is observed. We conclude that 
the source is accreting in the sub-critical regime. We discuss possible explanations for the double-peaked 
outburst. 

Subject headings: pulsars: individual (A 0535 +26) — stars: magnetic field — X-rays: binaries — X-rays: stars 



1, INTRODUCTION 

The Be/X-ray binary A 535+26 was discovere d during a 
giant outburst with ArieJ V (Rosenberg et al.ll 1 975b . At a dis- 
tance of d ~ 2kpc, it consists of the BO Hie optical com- 
panion HDE 245770 (ISteele et al.lll998h and a pulsating neu- 
tron star of spin period P sp i„ ~ 103 s, in an eccentric or- 
bit (e = 0.47) of period of P or b = 11 1.1 days (Finger et al. 
2006). The system presents different luminosity states asso- 
ciated with the activity of the Be star: quiescence, with X-ray 
luminosities L x < 10 36 ergs _1 (iRofhschild et al.ll2012h . nor- 
mal (type I) outbursts, with luminosities Lx ~ 10 36 ~ 37 erg s _1 , 
and giant (type II) outbu rsts, that can reach luminositie s 
Lx > 10 37 ergs~' (see, e.g. jFinger. Wilson & Harmonl 1996b . 
A 0535+26 presents two cyclotron resonance scattering fea- 
tures (CRSFs or cyclotron lines) in its X-ray spectrum at 
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-46 and -100 keV (see, e.g.. iCaballero et alj|2007T and ref- 
erences therein). Cyclotron lines are caused by resonant scat- 
tering of photons off electrons in the quantized Landau levels, 
and the cyclotron line energy is proportional to the magnetic 
field strength, £ cvc = heB/m e c = 11.6keV • B/10 12 G (e.g., 
ISchonherret"aT]|2007T) . 

Some accreting X-ray pulsars present a negative cor- 
relation between the cyclotron line centroid energy and 
the X-ray luminosity (e.g., 4U 01 15+643, V 0332+53, see 
Nakaii maetal.ll2006t iTsygankov eTaIll2006t iMowlavi et all 
2006). Other sources, like HerX-1, G X 304-1, and 
Swift J 1626.6-5 156 show the opposite trend (IStaubert et all 
I2007L Klochko vltini20ll IPeCesaretaU 120121) . As dis- 
cussed by iBecker et al.l (120121) . this bimodalitv is probably 
due to two different accretion regimes depending on the crit- 
ical luminosity L cr j t . For sources with Lx > L C nt (supercriti- 
cal sources), a radiative shock decelerates the infalling matter, 
with the height of the emission region increasing with increas- 
ing mass accretion rate, explaining the negative correlation 
between the cyclotron line energy and the X-ray luminosity. 
On the other hand, for sources with Lx < i C nt (subcritical 
sources), Coulomb interactions stop the infalling matter, and 
the height of the emission region decreases with increasing ac- 
cretion rate, e xplaining the positiv e correla tion observed (fur- 
ther details inlBecker et al . 1 120 1 2l see also iBasko & Sunvaevl 
119761 Staubert et al. 2007). Contrary to these sources, for 
A 0535+26 no significant correlation between the cyclotron 
line energy and the X-ray luminosity has been observed 
dTerada et al.ll20 06: Caballero et al. 2007), with the exception 
of a flaring episode during the onset of a normal outburst ob- 
served in 2005 during w hich the cyclotron line energy sig- 
nifica ntly increased (see Cab allero et al.ll2008t iPostnov et al.l 
2008 for details). These re sults refer to pulse phas e averaged 
spectroscopy; see however Kloch kov et al.l (1201 ll) for pulse- 
to-pulse analysis in which a positive correlation between the 
flux and the cyclotron line energy was observed. 

A 0535+26 showed a giant outburst and a series of nor- 

15 See, however, MulleretaT] j2012bl) . where it is shown that for 
4U01 15+64 the presence of a correlation depends on the model for the con- 
tinuum. 
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Fig. 1.— (a) Swift/BAT orbital light curve (15-50 keV) of the A 0535+26 
double-peaked outburst. The tick marks indicate the times of the INTE- 
GRAL, RXTE, and Suzaku observations presented in this work. The vertical 
dotted line shows the time of periastron. (b) Pulse period evolution derived 
from the RXTE observations. 

mal ones in 2009, all of them asso ciated with the peri - 
astron passage of the neutron star (jCaballero et al. 2009). 
The giant outburst reached a flux of F(i5-50)keV = (1.17 + 
0.04) counts s" 1 cm" 2 (~ 5.3Crat[3) at MJD 55179.0 and 
lasted about 40 days. In this work we focus on a double- 
peaked outburst that preceded the giant one. We present 
INTEGRAL, RXTE, and Suzaku observations of the out- 
burst. The Swift/BAT light curve of the outburst and the 
times of the observations are shown in Fig. Q] (a). The first 
peak of the outburst reached a flux of F(i5-50)k e v = (0.100 + 
0.004) counts s" 1 cirT 2 (-455 mCrab 16 ) at MJD 55048.0, that 
decreased to F(i5_50)kev = (0.047 ± 0.002) counts s" 1 cm -2 
(~214mCrab) at MJD 55051.0 and rose again, reaching 
^05-50)kev = (0.104 ± 0.005) counts s" 1 cm" 2 (~471mCrab) 
at MJD 55058.0 around periastron. The observations are de- 
scribed in Sec. [2] The results of the timing and spectral anal- 
yses are presented in Sec. [3] and the results are discussed in 
Sec.ffl 

2. OBSERVATIONS AND DATA ANALYSIS 

We made use of INTEGRAL data coll ected with the 
gamm a ray spectrometer SPI (20 keV-8 MeV, IVedrenne et al. 



2003), the imager IBIS (15 keV-10 MeV. lUbertirTeta! 



2003), and the X-ray monitor JEM-X (3-35 keV Xund et al 



2003). Three pointed observations were performed around the 
second peak of the outburst (MJD start 55054.89, 55057.99, 
and 55060.99). The JEM-X and SPI data analyses were per- 
formed using the standard analysis package OSA v9. For 
IBIS, we used OSA vlO, that contains a new energy calibra- 
tion for IS GRI 

RXTE dBradt. Rothschild & Swankll 19931) performed regu- 
lar pointed observations of A 0535 +26 during the outburst be- 
tween MJD 55046.86-55078.00 (Proposal ID P94323). We 
used d ata from the Prop ortional Counter Array PCA (2- 
60keV. |jahoda et aTT l996) and the H igh Energy X-ray Timin 
Experiment HEXTE (20-200 keV, Rothschild eTaD [jUT 
The RXTE data were analyzed using HEASOFT v6.7. We 
restricted the PCA spectral analysis to energies above 5 keV 
due to a feature around 5keV caused by instrumental Xe L 
edges dRothschild et alj|2006l) . 

Suzaku performed one observation of A 0535+26 dur- 

16 Fluxes derived from the daily Swift/BAT transient monitor results pro- 
vided by the Swift/BAT team. ~0.22 counts s cm -2 correspond to 1 Crab. 
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Fig. 2. — RXTE PCA and HEXTE energy dependent pulse profiles for 
one observation around the first peak (left, MJD 55047.86, ObsID 94323- 
02-03-01) and one observation during the decay of the outburst (right, MJD 
55065.84, ObsID 94323-03-02-01). The energy ranges are indicated in the 
left panels, and are the same for both observations. The PCA count rate is 
given in counts s" 1 PCU"'. Two pulse cycles are shown for clarity. 

ing the decay of the outburst (MJD start 55067.96, ObsID 
404054010). We used data from the t wo main instruments , 
the X-ray Imaging Spectrometer (XIS, Ko vama et al J 120071) 
and the Hard X-ray Detector (HXD, [Takahashi et all 120071) . 
The data analysis was performed using HEASOFT v6. 12 and 
CALDB versions 201 10913 for HXD, 20120209 for XIS, and 
20110630 for XRT We used XIS data between 0.5-10 keV 
for XIS 1, and data be tween 0.7-10 keV f or XIS 0, 3. XIS 
data were binned as in iNowak et al.l (1201 lb . Because of the 
XIS calibration uncertainties, especially around the Si K edge 
at ~1.8keV, a systematic error of 1% was assumed for the 
XIS data. The spectral analysis was performed with XSPEC 
V12.7.0. 

3. RESULTS 

3.1. Timing analysis 

We studied the pulse period evolution us- 
ing the phase-connection tech nique described by 



IStaubert. Klochkov & Wilms! (120091) . We extracted PCA 
light curves with 15.6 ms resolution in the 3.3^+2.8 keV 
range. Bary centric and orbital corrections were ap plied 
to the light curves, using the ephemeris from Fing er et al.l 
(2006). The evolution of the pulse period is shown in 
Fig. Q] (b). A decrease in the pulse period along the out- 
burst is observed. A linear fit can be used to describe the 
pulse period evolution between MJD 55046.89-55057.88, 
during the two main peaks of the outburst. We find a 
pulse period of P = 103.51263 + 0.00004 s and spin-up of 
P = (-0.459 + 0.005) x 10" 8 s s" 1 at MJD 55047.4. The pulse 
period during the decay of the outburst (MJD > 55057.88) 
can also be described with a linear fit. We obtain a spin 
period of P = 103.50886 ± 0.00001 s for MJD 55056.87 and 
a spin-up of P = (-1.104 ± 0.002) x 10" 8 s s" 1 . 

Making use of these pulse period values to fold the light 
curves, we studied the evolution of the energy dependent 
pulse profiles using all the RXTE pointed observations. Dur- 
ing the two main peaks of the outburst, the pulse profiles are 
remarkably stable, showing a complex structure at low ener- 
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Fig. 3. — Left; INTEGRAL spectrum of A 0535+26 (MJD 55057.99). (a) Data and best fit model, (b) Residuals of a fit including no absorption lines in the 
model, (c) Residuals of a fit including one line at ~ 46keV in the model. The filled circles represent the ISGRI data. The average flux during the observation 
^20-l0OkeV ~ 5.0 X 10 -9 erg cm~ 2 s _I . Right: Suzaku spectrum of A 0535+26 (MJD 55067.96). (a), (b), (c) have the same meaning as in the left panel. The 
black, red, and green symbols represent XIS 0, 1, and 3 data respectively. The average flux during the observation F20-l00keV ~ 0-5 x 10~ 9 erg cirT 2 s -1 , 



gies, and becoming simpler at higher energies. The shape of 
the pulse profiles at low energies slightly varies towards the 
decay of the outburst. As an example, the pulse profiles for 
one observation during the first peak of the outburst and one 
observation during the decay are shown in Fig. |2] 

3.2. Spectral analysis: 

We studied the broad band spectral continuum and its 
evolution during the outburst with RXTE, INTEGRAL, and 
Suzaku. It can be described by a powerlaw with an expo- 
nential cutoff, typical for ac creting X-ray binaries (see, e.g., 
IWhite. Swank & Holtl [19831) . We used the XSPEC model 
cutoffpl, given by F(E) oc E~ T e~ EIEmi , where F and E Md 
are the photon index and folding energy respectively. 

The Suzaku spectrum was modeled including an additional 
black body component, photoelectric absorption and a Gaus- 
sian line in emission to account for the Fe Ka line. The 
presence of the Fe Ka line is marginal, and the line en- 
ergy and width have been fixed to 6.4 and 0.5 keV respec- 
tively. Its inclusion in the model yields an improvement in the 
*L/d.o.f from 1.43/581 to 1.37/580. The equivalent width 
of the line is EW (Fe Ka')=70 + 2 eV. The inclusion of the 
black body component in the model leads to an improvement 
in the x^/d.o.f. from 3.01/582 to 1.37/580, giving an F-test 
significance > 99.99%. The best fit value of the black body 
temperature is ksT — 1 .26^ ^ ke V. Setting the cr oss-sections 
and the abundances to those from Ve rner et al.l (1 1996b and 
IWilms. Allen & McC ray ( 2000) respectively, we obtained an 
absorption column of Nn = 0.70 + 0.03 x 10 22 atoms cm -2 . 

To model the RXTE spectra, we added a Gaussian emis- 
sion line to account for residuals present in the 5-7 ke V range. 
The energy of the line was left free, and the width of the line 
was frozen to the best fit value in each observation. The line 
energy obtained for the different observations varies between 
5.0 and 6.6 keV, consistent with the instrumental Xe L edge at 
~ 5 keV and the 6.4 keV Fe Ka line. It was not possible to 
model the two components independently. 

In the INTEGRAL, Suzaku, and most of the RXTE ob- 
servations, a significant absorption-like feature is present in 
the residuals at E ~ 46keV, that we model using a Gaus- 
sian optical depth profile t(E) = Te^ E ^ E ^' ^ 2o ~ ) , which mod- 
ifies the continuum as F'(E) = F(E)e~ T( - E> . As an example, 
the broad band spectrum of the second INTEGRAL observa- 



tion is shown in Fig. [3] (left). In this case, the inclusion of 
a cyclotron line in the model yields an improvement in the 
xfjd.o.f. from 8.8/132 to 0.88/129. The Suzaku broad band 
spectrum is shown in Fig. [3] (right). The inclusion of a cy- 
clotron line in this case improves the^ 2 d /d.o.f. from 1.73/583 
to 1.37/580. 

We studied the evolution of the cyclotron line centroid en- 
ergy, photon index, and folding energy as a function of time 
and X-ray luminosity. To study the evolution of the different 
spectral parameters, we selected the first 16 RXTE observa- 
tions (MJD 55046.86-55065.84), which are the ones where 
the cyclotron line was significantly detected. The first har- 
monic line at E ~ lOOkeV, observed in brighter outbursts 
of the source (e.g.. lCaballero et aLll2007l) . is not significantly 
detected in the observations presented here. The results are 
shown in Fig. [4] No dependency of the cyclotron line en- 
ergy with the X-ray luminosity is found. We checked that the 
choice of continuum model does not affect the line param- 
eters by fitting the continuum using an expo nentially cutoff 
power law, a Fermi Dirac cutoff power law (Tanaka 1986), 
and a power law with an exponential cutoff starting at a cutoff 
energy (the highecut model of XSPEC). In all three cases 
the line parameters rema in the same within their error bars 
(see also Caballero 2009). A linear fit to the E cyc vs. Lx val- 
ues gives a slope of (-3 + 3) keV/10 37 erg s _I using RXTE 
data only, and (-0.9 + 2.0)keV/10 37 erg s" 1 using RXTE, IN- 
TEGRAL, and Suzaku data. We included inter-calibration 
constants for the different instruments in the separate RXTE, 
INTEGRAL, and Suzaku spectral analysis. For the Suzaku 
analysis, we fixed the constant to 1 for XIS and allowed 
the others to vary. In the case of the GSO detector of the 
HXD, we fixed the constant to cqso = 0.88 (the value ob- 
tained before adding the cyclotron line to the model), because 
when unconstrained, the best fit results in an unrealistic value 
of cqso - 0.4 + 0.1. By doing this, the folding energy de- 
creases from Efou = 44*' to 31 + 2keV, while the pow- 
erlaw index and the cyclotron line energy remain constant 
within the errors. As seen in the middle and bottom pan- 
els of Fig. [4] the photon index significantly decreases with 
increasing luminosity. The folding energy remains constant 
during the main part of the outburst, down to a luminosity of 
£(3-50) keV ~ 0.26 x 10 37 ergs~', and shows an increase in the 
last observations of the outburst. Both the folding energy and 
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Fig. 4. — Left: Evolution of the cyclotron line centroid energy £ C y C (a), 
photon index T (b), and folding energy £f i,i (c) with time. In panel (a), the 
crosses represent the Swift-BAT light curve in arbitrary units. Right: evolu- 
tion of £ cyc (d), r (e), and £f„id (f) with the (3-50) keV X-ray luminosity. The 
RXTE, INTEGRAL, and Suzaku observations are represented with triangles, 
squares, and a circle respectively. In panel (d), the dashed line represents a 
linear fit to the data. Error bars are at 90% confidence. 

the photon index increase with decreasing luminosities. The 
values plotted in Fig.|4]have been obtained leaving these two 
parameters free in the spectral fits. The increase of the pho- 
ton index with decreasing X-ray luminosity indicates a soften- 
ing of the spectrum, while the increase of the folding energy 
with decreasing X-ray luminosity indicates a hardening of the 
spectrum. In order to check if the spectrum becomes softer 
or harder as the luminosity decreases, in the last three RXTE 
observations we fixed the folding energy to the mean value of 
the first 13 observations, Ef a n = 18.4keV. By doing this, a 
slope of (-0.44 + 0.04)/10 37 ergs~' is obtained from a linear 
fit to the T vs. Lx values, with a Pearson correlation coeffi- 
cient of 0.82. This result indicates a softening of the spectrum 
with decreasing luminosity. 

4. DISCUSSION 

We have presented the first observations of a double- 
peaked outburst of A 0535+26 with RXTE, INTEGRAL, 
and Suzaku. The double-peaked outburst was the precur- 
sor of a giant one in December 2009. Two double-peaked 
outbursts were also observed b efore a giant one in 1994 
(Finge r/Wilson & Harmonl 19961) . w hile other giant ou tbursts 
do not show any precursor outbursts dTueller et al.l2005h . An- 
other outburst with a peculiar shape (two separate peaks) was 
observed two orbital phases after the 2009 giant outburst (see 
iCamero-Arranz et al.ll2012l) . 

We measured a pulse period of P = 103.51263 + 0.00004 s 
and spin-up of P = (-0.459 + 0.005) x lO^ss" 1 at MJD 
55047.4, and a spin-up of P = (-1.104 + 0.002) x 10~ 8 s s" 1 
at MJD 55056.87. These results are i n agreement with the 
values obtained with Fermi-GBM by Cam ero-Arranz etaD 
(120121) . who also present the long-term pulse period history of 
A 0535+26. We note that the changes in the spin-up rate along 
the outburst could be due to the outdated orbital solution. A 
new orbital solution is required to explore the pulse period 
evolution in more detail. A spin-u p during normal ou tbursts 
has also been observed in the past (Caba llero et al.l2008l) . pro- 
viding evidence for an accretion disk around the neutron star 
during normal outbursts. The energy and luminosity depen- 
dent pulse profiles are remarkably stab le compared to past 
observations of the source (see, e.g., ICaballero et al.l 120071 



and N aik et al.l2008l). and consistent with tho se obtained with 
Fermi-GBM bv lCamero-Arranz et al.l (120121) during the same 
period. 

We observe no significant variation of the cyclotron line en- 
ergy with the X-ray luminosity during the double-peaked out- 
burst, and a significant decrease of the photon index F (spec- 
tral hardening) with increasing luminosity. A similar behav- 
ior has been observed in other outbursts of the source (e.g., 
ICaballero et~aiT2 008; Mu ller et al.l2012al) . In accreting X-ray 
pulsars, a positive or negative correlation of E cyc with Lx is re- 
lated to different accretion regi mes. Following the theoretical 
predictions from iBecker et al.l (120121) . the critical luminosity 
for A 0535+26 corresponds to L crit ~ 6.78 x 10 37 erg s" 1 . The 
observations presented here are well below that limit, with 
i(3-50)kev ranging between ~ 0.04 and ~ 0.56 x 10 37 ergs~', 
and therefore A 0535+26 is probably accreting in the sub- 
critical regime during the outburst. A positive correlation be- 
tween E cyc and Lx is expected for sub-critical sources. The 
fact that no correlation is observed might be due to the fact 
that our observations are near or below the Coulomb stopping 
limi t, L m ,i], where little variation of E cyc with Lx is expected 
(see|Becker et al. 2012, for further details). 

The double-peaked shape of the light curve could be due 
to p erturbations in the Be disk around the optical compan- 
ion. iMoritani etail (1201 lb showed that A 0535+26 exhibited 
strong Ha variability during, before, and after the giant out- 
burst in 2009, and suggested that strong perturbations in the 
Be disk started about one cycle before the giant outburst. 
In addition, the Ha EW and V-band brightn ess showed an 
anti-correlation before the 2009 g iant outburst (lYan. Li & Liul 
l2012tlCamero-Arranz et al.ll2012l) . This indicates that a mass 
ejection event took place before the giant outburst, pro- 
ducing a low-density region in the inner pa rt of the disk 
that c ould explain th e double-peaked profi le (lYan. Li & Liul 
120121) . Interestingly, lYan. Li & Liul (120121) reported a simi- 
lar Ha EW and V-band brightness evolution before the 1994 
outburst, when tw o double-peaked outbursts took place be- 
fore the giant one (Finge r. Wilson & Har mon 1996). Such an 
anticorrelation was not observed before the 2005 giant out- 
burst, that took place without precursor outbursts. Double - 
peaked light curves have recently bee n observed in other 
sources, for instanc e GX 304-1 (jNakaiima et al.l 120121) and 
XTE J 1946+274 dMiiller et al.1 120121) . Whilethe spectral 
shape again remained relatively constant over the outburst for 
the latter, this is a special case, since two outbursts per or- 
bit are observed. In the case of GX 304-1 the double-peaked 
outburst took place before a giant one, similar to what has 
been observed for A 0535+26, suggesting that double-peaked 
outbursts could be i ndicators of upcoming giant outbursts 
(Nakaiima et alj 120121) . The double-peaked outburst in the 
case of GX 304-1 was similar to the one of A 0535+26 in 
terms of duration, lasting about 40 days. While in the case of 
A 0535+26 the two peaks reached a similar flux level, in the 
case of GX 304- 1 the second peak reached a flux two times 
higher than the first one. The intensities of the two peaks were 
^05-50)kev = (0.098 ±0.004) counts s _1 cirT 2 (~445mCrab) at 
MJD 56075 andF (15 _50)keV = (0. 196 + 0.008) counts s" 1 cirT 2 
(~891mCrab) at MJD 56087, with a flux between the 
two peaks that dropped down to f(i5_50)kev = (0.021 ± 
0.001)countss-' cm' 2 (~95mCrab) at MJD 56082. Note 
that the distances of the two so urces are rather similar, 
with G X 304-1 being at ~ 2.4 kpc dParkes. Murdin & Mason 
1980). The giant outburst of GX 304-1 lasted about 50 days, 
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showing a first weaker peak that reached a X-ray flux of 
^(i5-50)kev = (0.163 + 0.007) counts s^ 1 crrT 2 (~741mCrab) 
at MJD 56208. The flux then dropped down to zero (for- 
mally F ( i5-50)kev = (0.004 + 0.021) counts s" 1 cnr 2 at MJD 
56223), and increased again, reaching F (i5_50)keV = (0.406 ± 
0.015) counts s^cnT 2 (~ 1.85 Crab) at MJD 56234. Fur- 
ther simultaneous X-ray and optical observations of double- 
peaked outbursts are needed to better understand their origin 
and relation to giant ones. 
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